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1) Calcul d’une valeur approchée de In(2) 


a) Approximation de In(2) par le produit de Seidel 


Théorème : formule du produit de Seidel ! 


e Pour tout réel x > 0, on considère la suite (u, (x)) définie pour tout n € N° par : 


neN* 


n 


_ RE 2 2 2 
HE Née 1x2 1x4 1e * Xi Eve 
k=1 


2 2 2 , 2 
1+vx VAE Er 
AR 


: In x 
Alors la suite (u,) converge vers : 


| In x 
M Un Q) = 


l _—_ 
(six = 1, on remplace — par sa limite 1). 


e Six = 2, la suite (u,(2)) définie pour tout n € N* par : 


neN* 
T D... 7 2 2 2 
ROSE tt 
k=1 
— 2 x 2 X - X … X : 
Ti su. TE oUeT 
1+VV2 1+.,/VV2 


est strictement décroissante et converge vers le logarithme népérien de 2 : 


lim Un = In2 & 0,6931471806 
n— +00 


Démonstration 
http://www.les-mathematiques.net/phorum/read.php?2,244589 
http://people.math.binghamton.edu/dikran/478/Ch6.pdf ($6.1.3) 


Soit x € R*. On introduire son cosinus hyperbolique et son sinus hyperbolique, définis par : 


" e*+e * 
COsh x = ——— 
2 
| e*—-e * 
sinh x = 
2 


1 Philipp Ludwig von Seidel (1821-1896), mathématicien et physicien bavarois, a montré cette formule en 1871. 
https://fr.wikipedia.org/wiki/Philipp_ Ludwig von Seidel 
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e2+e 
x 24% E 2% E 2% _ 2 
x sinh> 
= cosh — X x 
2 
sinh x x Sinh 7 
= cosh 5 X x 
2 
En remplaçant x par = , on obtient : 
+ . 1 À 
sinh 5 x Sinh 72 
Æ cosh=; X FA 
2 22 
En remplaçant x par : , on obtient : 
À | 
sinh 72 x  Sinh 53 
: = cosh 33 X FA 
En remplaçant x par = , on obtient : 
sinh sinh 
23 X 24 
 — cosh 74 A 
p. 2? 
Plus généralement, pour tout n € N, 
: s.* , x 
sinh Sn x sinh SR 
X = cosh pn+1 X X 
on pn+1 
En effet : 
x x 2X 2x 2 2 
X e27t — e 27 e2nri — e 21 XX __X 
mé ner GT (in) (en) 
ÉLUS 2 _ 2 - ( 
X X n 2% - 22% 
PAL PAL pn+i pn+i 
X 
x x x x onFI _ po onFi 
(er JE er) (er _ er) RE n QT e2 e 2 
_ 2 X 2x CT 
pn+1 pn+1 
x 
X sinh sr 
= cosh Jai Ne — — 
pn+1 
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x x sinh 7 
= cosh — X x = cosh> x cosh>3 x > — 
57 
x x x sinh 33 
cosh — X cosh — X cosh — x ——— 
2 22 23 +. 
23 
—— 
X X x X sinh 57 
cosh — X cosh — X cosh — X cosh — X —— = 
2 2? 23 24 X_ 
24 


X n 


sinh DT x 
Sn: [| cos 


pn k=1 


Pour tout n € N, 


1 x sinh _. 


+ 
x sinh TA 


& 
Ï 
| 
in 
Éa 


lim | — 
n+® | sinh x = n-+00 cosh -— 


x li [| 1 
, = lim — 
sinhx n-+0 cosh 37 


k=1 


Pour tout réel x > 0, il existe un unique réel 8 tel que x = n8 & e* = 6. 
6 D-0 0-1 
2 2 20 


sinh x = 
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X In 0 _20Mm0 28 In60 
sinhx 2-1 2-1 (98—-1)(8+1) 


20 
en. FA RS LT à 
h x e*+e 2 (e*)*X+(e-*)2%X 92 +(9"1)2% gr +09 2 
D pi. VU UT e uDT  GE 
1 1 1 1 1,17 1271 2 1 
002% +02%0 2% OX Æ+0X 2% GA+O pK 1+1 
. 1 = Es _ LE Le 
202* 202 202k 202% 
JL 
1 202* 


On a vu que 


de 


20 In 8 , II 202% : il x] | 2 ) 
HDŒED ve eE. | jo CL 
G-D@+D née(ll ne, nel 


n 

1 1 1 1 1 1, 1 1 1 n 1 
| | O2 = O7 x 02 x O2 x …. x O2 = Qui artrst te = GRR 
k=1 


. : : ; ris . 1 
On utilise la formule donnant la somme de n premiers termes d’une suite géométrique de raison ; 


1 
: 1 
et de premier terme Ô =, 


D DD termes D () 


2 


26 In 8 [I 2 
————————__ — = X ——————— 
D=DOED 44 _ 


k=1 O2k 1 +71 
5. à = 2 . 2 
= jm, (8) « in 6 x im (T1) 
n—+00 n—+00 n—+00 
k=1 021 +7 k=1 O02K-1 +1 
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Dans le produit | 


k=1 O2k-1 + 1 


Lorsque k varie de 1 à n, p = k — 1 varie de Oàn — 1. 


n-1 — 
2 2 2 2 
Ts | NE 
k=1 02-141  p=0 er 1 02+1p-102P +1 
he 
ll. 
= | | —— 
ne p=1 02P +1 
-1 
26 In 8 —. 2 I 2 ox? 
sn = im | —— = 
(8 — 1)(8 +1) n=+o 0+11l6%41 0 +1 
n—-1 n 
In 8 h [| À h [| 2 
— E Le 1 - aus 1 
p=1 1 + 02P Do p=1 1 + 02P 
ee est la limite quand n — + du produit : 
n 
2 2 2 2 


1 
p=1 1 + 02P 
2 2 2 


cn (Te 
n 1 


, on effectue le changement de variable p = k — 1. 


n—-1 


p=1 02P +1 
+00 
2 
1 
p=1 1 + 02P 
2 


Tire *i+e "ire" *LE er 


2 


EX À  ———Ù KE 
14% = T + 672 
L+VV8 7 


() 
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b) Approximation de In(2) par les sommes de Riemann des 
1/(n+k) pour k allant de0àn-1oudelàn 


Définition 

Soit a et b deux réels tels que a < b. 

Soit f une fonction continue par morceaux sur [a ; b] à valeurs dans R. 

Pour tout n € N°, on appelle sommes de Riemann d’ordre n associée à f les sommes : 


n-1 


buts 


LL Lie. 


En 7) 
{a REP EC EE Ne +f(atn—)] 


Théorème : convergence des sommes de Riemann 
Les suites (r }nen” et (S)hnen* convergent vers l’intégrale de f entre a et b. 
b 


li % = Jim Se J roa 
nn +00 


Démonstration 


Voir p.635 dans MPSI, mathématiques tout-en-un, Dunod, 4° édition, 2015 


nn ue b : 

Théorème : approximation de In a par des sommes de Riemman 

Soit a et b deux réels tels que 0 < a < b. 

On considère les suites (r, }nen* €t (S)nen* définies pour tout n € N* par: 


RS 1 1 1 
LR D rues ur A tri 
a a 
CNT O1 1 1 1 1 
Fa +k p—+1 a 2 +3" ju: +n 
a a À 47 a 4 


| _ b 
Alors les suites (7, )nen* et (S)hen* convergent vers le logarithme népérien de — 


b 
lim 7, = Jim Sn =In-— 
a 


n—+oo 
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Démonstration 


Soit a et b deux réels tels que 0 < a < b. 
2. Pa 1 4 # . 1 
On applique le théorème ci-dessus à la fonction inverse f : x > f(x) = T 


La fonction inverse f est continue par morceaux sur [a ; b] à valeurs dans R. 


Pour tout n € N°, les sommes de Riemann d’ordre n associée à f sont les sommes : 


n-1 n-1 n-1 b—-a n-1 
| +k2—) D) 1 nn 1 
M = a = —— — ————————— = —————— 
b—-a ji b—a > n 
a — jé É K=oatk—— Koatk— 04p-atÀ 
n-1 
: 1 : 1 . 1 . 1 1 
Li —tk —+1 7 —+2 —+n-1 
_— 2 | dl 
a a a a 
b — > ne b — ) > 1 
RE — = 
k=1 k=1h 
| 
a 
: 1 1 ” 1 Fr 1 
+1 2 +3 p—+n 
En | __i 41 
a a a a 
Les suites (7, )nen” et (S)nen* convergent vers le logarithme l’intégrale de f entre a et b. 
b b 
1 b 
lim r, = | roc _ ET = [ntJ}=Inb-Ina=In- 
n—+00 t a 
a a 
Théorème : approximation de In x par des sommes de Riemman 
Pour tout réel x > 1, on considère les suites (rm Ce et (sa Ho) 
définies pour tout n € N° par : 
n-1 
1 1 1 1 1 
A D —— © 5 — + 5 — + 5 — ++ 
K=0 qui r—1" x ti >=" +2 >= tn-i 
n 
1 1 1 1 1 
Sn(x) — D ——- PR Cp à + 
ki =qu+tk =jn+i "+2 >=" +3 juin 
Les suites (rm ce et (sa œ)) convergent vers le logarithme népérien de x. 
lim r,(x) = lim s,(x) =Inx 
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Démonstration 
Soit un réel x > 1. On applique le théorème ci-dessus avec a = 1 et b = x. 


Les réels a = 1 et b = x vérifient a < b puisque 1 < x. 


On considère les suites (7 )nen* €t (Sn)nen* définies pour tout n € N* par : 


n-1 n-1 n-1 


Go) 5 1 1 1 
= LL 
a 1 
EE : 
nn ñ 7] ñ 
x—1 zx=-1"1 >x-1"2 x stat 
= 1 = 1 
Sn (x) = 7] = 
k=1 k=1 x tk 
a 
EE 
nn En En 
RE SE | x=1t2 x=1t3 LL 


Les suites (r, 6) 


: _ ss 
Jim mix) lim sr) 1 — In x 


Exemple si x = 2 


Les suites (r,(2)) et(s, ce. définies pour tout n € N* par : 


neN* 


n-1 n-1 n-1 


D 


k=o5=qu+tk koqntk ko 
: 1 1 1 __— 1 
[on n+1 n+2 n+3 n+tn—-1 


n 


@=Y 1 : 1 F 1 1 n 1 
Sn An+tk n+1 n+2 n+3 n+n 


convergent vers le logarithme népérien de 2. 


lim r,(2) = lim r,(2) =In2 = 0,6931471806 
n—+00 n—+00 
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Démonstration de l’exemple si x = 2 sans utilisant le théorème ci-dessus 


https://www.ilemaths.net/sujet-une-suite-qui-converge-vers-In2-119013.html 


On considère la suite (u, }nen* définie pour tout n € N* par : 


_ 


ÿ :  . 
1: = — … 
" Antk n+1 n+2 n+3 n+n 


Calculer u; ; U2 ; U3. 


a) Démontrer que pour tout n € N*, 
1 


Mat Gr DCn+D 
b) En déduire les variations de la suite (u,). 
On considère les fonctions f et g définies sur ]0 ; +oo[ par : 


1 
fG)=mx-(1-=) et g(x) =Inx -(x-1) 
a) Etudier les variations de f et g. 


b) Démontrer que pour tout réel x > 0, 1 — = Xe xl 
pt : | di 4 
c) En déduire que pour tout p € N°, — < In RE 
p+1 D p 


a) Soit n € N°. Écrire l’encadrement de la question 3.c) en donnant successivement à p 
les valeursn,n+1,n+2,...,n+n. 
b) En sommant l’encadrement de la question 3.c) pour p allant de n à 2n — 1, montrer que 


1 
Un SM2£Suts 


Démontrer que la suite (u,) converge vers In 2. 
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Solution 


1. 
1 
D L 1 
AT ask ni 
k=1 
sS 11 , 1 _n+t2tnti  2n+3 
27 Lntk n+1 n+2 m+Dh+2) G+DM+2 
3 
ss En + À 
8 Lintk n+l n+2 n+3 


_(+2)(n+3)+(n+1)(n+3)+{(n+1)(n+2) 
Es (n + 1)(n +2)(n +3) 


_M PORT G ER PARIS tin E2 3n2 + 12n +11 
E (n + 1)(n + 2)(n + 3) _(n+1)(n +2)(n +3) 
2. a) Pour tout n € N*, 
n+1 
_ à 1 
Sn PE TES nn 


Sur la 1°® somme, on effectue le changement de variable : p = k + 1. 
Lorsque k varie de 1àn+1,p=Kk+1variede2àn+2 


n+2 n 


1 1 
Uni Un = D _. 
p=2 P k=1 
n n 
D 1 . 1 1 _. 
E int+p ntn+i ntez _n+1l : 
p= =? 
n n 
> 1 = 1 
Zintp aintk 
1 . 1 1 1 : 1 2 
UM nl 2n+2 NE 2nt1 2+1) 2R+1 
1 1 _2RRDSOREI 2m4+27-2n=1 


2n+1 2h+1  2n+D@n+1) 2n+1)@2n+1 
1 


"2n+DQ@n+1) 
Pour tout n € N°, 
1 
Ut Gr DCn+ED 
b) Pour tout n € N*, 
1 
Ha M EDS D 
Un+1 Un > 0 >  Unsi > Un 


La suite (u,) est strictement croissante. 
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3. a) Les fonctions f et g définies sur [0 ; +oco[ par : 


ft = mx-(1-") et g(x) =Inx -(x-1) 


sont dérivables sur ]0 ; +ool. 


ei 1 x 1 _x—1 aus &t 1) 
f(x 7,272 42 2 (stdusigne de (x 
1 1 x 1—-x | 
g'x)=--1=---- est du signe de (1 — x) 
* x + # 


f@=mi-(1-;)=0-0-0 
g{i)=M1-(1-1)=0-0=0 


x 0 1 + oo 
signe de f'(x) — + 
+00 0 +00 
variation de f LR 
0 
x 0 1 + 00 
signe de g'(x) + 0 _ 
0 
variation de g D 
— OO — O0 


nr 


b) Pour tout x E [0 ; +oo!, 
1 1 
fG)=mx-(1-=)20 = Inx>1-- 
x x 
gx)=Inx-(x-1)<0 = Inx<x-1 


1 
1—--<Inx<x-1 
F4 
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——. 1 
c) Pour tout p € N six == 1+5,alors 


1 
RS 
P 


et 
=. 
x p+1 p+li p+li p+l 
L’encadrement 


1 
1—--<Inx<x-1 
4 


devient 


< — 
br p p 


4, a) Soit n € N°. Écrivons l’encadrement de la question 3.c) en donnant successivement à p les 
valeursn,n+1,n+2,...n+n—1. 


1 ,m+t1,1 
sip=n, < In < — 
P n +1 n n 
1 n+2 1 
UE HE ad nil 
 _ L > mi 7 
ne DES Aad2 nie 
1 2n 1 
sip=n+tn-1=2n-1, — < In 


b) Si on somme les encadrements 


1 1 1 


p+1 p p 
pour p allant de n à 2n — 1, on obtient : 
2n—-1 


2n—-1 2 
1 p+l1 1 
Dr Dm 
mp+li _ p p 
p=n p=n 


e Dans la somme de gauche on = , on effectue le changement de variable k = p + 1 —-n de 


sorte que p + 1 = n + k. Lorsque p varie den à2n—-1,k=p+1-nvariedn+1-n=1 
à 
2n—1+1-1n="n, 

2n-1 


1 UE — + + ++ 
p+i Zntk  n+I n+2 n+3 n+n 
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- 
Fe mt à jé 
p=n 


n - 


e Dans la somme de droite , on effectue le changement de variable k = p — n de sorte 


A 
= 


2n-1 


DE >  — BE  d 2. 
n+k n+0 n+tn ZLintk on Re D nu 


e La somme centrale est : 
2n-1 2n-1 2n-1 2n-1 


Dis = me +D- In p] = 2H 


Dans la somme >%2,! In(p + 1) , on effectue le changement de variable k = p + 1. 


Lorsque p varie de n à 2n — 1, k = p + 1 varie de n + 1 à 2n. 


2n—-1 2n 2n—-1 2n-1 2n—1 
p + 
> Mn > In k — Due Fa In k + In(2n) —- Inn -— Inp 
p=n k=n+1 k=n+1 p=n+1 
2n-1 2n—=1 
> In k = > np 
k=n+1 p=n+1 
2n-1 +1 
> In? = In(2n) — Pre In2 
p n 
p=n 
2n-1 
> hé 
p=n 
L’encadrement 
2n—1 2n—1 2n-1 
1 p+1 1 
| 
mp+1 Z p — D 
p=n p=n p=n 


devient 
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5. Pour tout n € N*',ona 


L 
Un SM2£uU+SE 


1 1 
< — > — — < 
In 2 Un ts , +m2<u, 


1 
on tM2<u < In2 


1 1 1 
lim —-=0— lim (-)=0— lim (-5+n2) = m2 
2 2n 


n—+0o n— +00 n n— +00 


lim In2=0 


n—+00 


La suite (u,) est encadrée par deux suites qui convergent In 2, donc d’après le théorème des 
gendarmes, la suite (u,) converge vers In 2. 


lim u, =In2 
n—+00 
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c) Approximation de In(2) par la somme de Brouncker de 


1/[(2k-1)(2k)] pour k allant de 1àn 


Théorème 
e Pour tout réel x € [0 ; 1], on considère la suite (by (x)) 


définie pour tout n € N° par : 


neN* 
7X 2k-1 2k 
X X 
b > 1Y + D (- 
n Cu Ces ) 
a 5 4 x? x° x27-1 x22 
à _— 


23 4 5 à ‘on L 2 
Alors la suite (b, (x)) converge vers le logarithme népérien de (1 + x) : 


lim b, (x) = In(1 + x) 
n—+00 


e Six = 1, la suite (b,(1)) définie pour tout n € N* par : 


neN* 
2n-1 | n 
b. (D = Dev 1 = 1 Li À 1 
F7 Li+1 AUk-1 2k/ 2 3 4 5 6 2n—1 2n 
< 1 1 1 1 1 


Ë GK DXZR 1x2 3x4 5x6" Ÿ(n-1)x2n 


est strictement croissante et converge vers le logarithme népérien de 2 : 
lim b, (x) = In2 = 0,6931471806 
n— +00 


Remarque : on montre en annexe p.40 que la suite b,, (1) est égale à la somme de Riemann 
n 


1 L L 1 
ne 


Zntk n+1i n+2 n+3 n+n 


Démonstration 
e 1° point 
Pour tout x € R \ {1}, on sait que la somme des (n + 1) premiers termes d’une suite (u, )nen 
définie par u, = x", qui est la suite géométrique de raison x et de premier terme 1, vaut 
Le 1—x7+1 1 xtt1 


k 2 3 n 
LL ER RE  —_—_—_—_—_—_————_ — — 
» 1—-x 1—x 1—-x 


k=0 
1 1 n+1 
=D xt+ (E) 
k=0 
x+1 
= XT — = xeE(x 
1—-x 1 — @) 


où la fonction &, définie sur R \ {1} par £(x) = —_ tend vers 0 lorsque x tend vers 0. 
lim (x) = 0 
X—0 


n 
1 
= D x +rre() =L+x txt + x ++ xt + ate(x) 
k=0 


: di . 1 
C’est le développement limité en 0 à l’ordre n de la fonction x — 2. 
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En remplaçant x par —x , la relation (E) devient, pour tout x € R \ {—-1}, 


n n 
1 —X n+1 —1 n+1,n+1 
= Do + (C2 — D CDKxt + (Co D M 
1+x 1+x 14% 

k=0 k=0 


—1 n+1,n+1 _—1 n+1, 
ee EE = 2000 


_syn+1 
où la fonction 6, définie sur R \ {—1} par 8(x) = = tend vers 0 lorsque x tend vers 0. 


lim 8(x) = 0 
X—0 


n 
1 
au — D CDKxt + x 0x) = 1-x+x7 x + + (1) + x (x) 
k=0 


’  . à . 1 
C’est le développement limité en 0 à l’ordre n de la fonction x — 


; +. : 1 
Le développement limité en 0 à l’ordre n — 1 de la fonction x Ts St: 


1 
1+x 


= D C-DExt + (Da + 2700) 
k=0 


(—1) xt + x0 (x) = x 1[(—1) x + xO(x)] = x Lo(x) 


où la fonction 6, définie sur R \ {—1} par o(x) = (—1)"x + x8 (x), tend vers 0 lorsque x tend 
vers 0. 


lim o(x) = 0 
x—0 


n-1 
1 
FT D CDixt LE To (0) = 12 +7 oO EC) 
k=0 


Sur ]—1 ; +oo! , la fonction x + In(1 + x), qui est une primitive de la fonction x — , admet 
en 0 un développement limité à l’ordre n qu’on obtient en prenant les primitives de chacun des 
termes du développement limité à l’ordre n — 1 de x — , auquel on ajoute la constante 
d’intégration In(1 + 0) = In(1) = 0. On obtient : 


n-1 


In(1+x) = ) (—-1)* 


+1 x2 x 4 


ar 2 7 4 Ca a) 
RER tr 23 4 es 


où la fonction t vérifie 
lim r(x) = 0 
x—0 


C’est le développement limité en 0 à l’ordre n de la fonction x + In(1 + x). 
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Le développement limité en 0 à l’ordre 2n de la fonction x + In(1 + x) est : 
2n-1 


In(1 + 0 = 20 FEES r + 28e) 


i+1 

x2 x x x27-1 x22 

= X——+——-—+e+(—-1)772 + CD + x MT (x 
273 4 COS SR FENTE T() 
x? x _x* x27-1 x27 , 

SR CE 

D A nel 2n 0 

x2k71 x2k 

— a 2n 

es )++ 14) 


Pour tout x E[—1:1[,x2 € ]0:1[, donc 
lim. x = lin (x°)7 = 0 
n—+00 n—+00 


= girl Le 2k-1 2k 


i _ ds | 
Fu 2. D ae) (= ) = na +0 


Ll= 


S1 on fait tendre x vers 1 par valeurs inférieurs, on obtient par continuité : 
2n—1 


n 
1 1 
i — : — 
Jim | (— D'- = im, GC -—)= In(1 + 1) = In(2) 


i=0 k= 


Conclusion : pour tout réel x € [0 ; 1] , on considère la suite (by Go). définie pour tout n € 


N° par : 
— 2k-1 2k 
X x 
— i = 
Ba) = 2 C o D (- x) 
x? "À 4 x? x° x2771 x27 
1728 2 = 6 ‘Mel 2 
Alors la suite (b,,(x)) converge vers le logarithme népérien de (1 + x) : 
lim b, (x) = In(1 + x) 
n—+00 
e 2° point 
Si x = 1, la suite CD) définie pour tout n € N° par : 
2n-1 . n 
RON 1 —)=1 pe — . F 1 1 
‘  Ai+1 AUxk-1 2kJ 2 3 4 5 6 2n—1 ?2n 
1= = 
1 1 1 1 1 


_ GK DXZ2R 1x2 3x4 5x6" Ÿ(n-1)x2n 


est strictement croissante et converge vers le logarithme népérien de 2 : 
lim b, (x) = In2 & 0,6931471806 
n— +00 


1 1 2k 21 2k=2k+1 1 
2k=—1 2k (2k—1)X2k (2k—1)Xx2k (2k—1)Xx2k (2k—1)x2k 
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Méthode de Brouncker 
La convergence de b, (1) vers In 2 peut se retrouver par la méthode de Brouncker. 


Pour tout réel x > 0, on considère l’intégrale 
1+x 


1 
I(x) = l dt = [Int}i** = In(1 + x) —In1 = In(1 + x) — 0 = In(1 + x) 
1 
En particulier, si x = 1, 


1(1) = fra = In(1 + 1) = In(2) 


L’algorithme de Brouncker’ consiste à calculer une valeur approchée de l’intégrale 1(1), vue 
comme l’aire de la surface délimitée d’une part par l’axe des abscisses et la courbe (hyperbole) 
représentative de la fonction inverse, et d’autre part par les droites verticales d’équations 
respectives y = 1 et y = 2. Cette aire se calcule approximativement en sommant les aires de n 
rectangles R;, R;, R3, ..., R, sous l’hyperbole. 


Pour tout k € [1 ; n] , chaque rectangle R, a pour longueur (ou base) L}, largeur (ou hauteur) £z, 
et aire Ay —= Le X L3. 


On construit les rectangles en respectant deux règles : 
2 : à 1 
e Chaque rectangle R, a son sommet supérieur droit S, sur l’hyperbole, donc S} (xx 3 Yr = =). 
k 
e Chaque rectangle R, a sa longueur L}, qui vaut la moitié de la longueur du rectangle sur lequel 
il est posé. 


Chaque étape de construction consiste à poser, sur chacun des rectangles construits à l’étape 
précédentes, deux rectangles adjacents de même longueur : un rectangle à gauche de largeur non 
nulle, et un rectangle à droite de largeur nulle. Aïnsi, pour tout k € [2 ; n], il y a à droite du 
rectangle R, de longueur L, et de largeur £, # 0 un rectangle R? de même longueur L, et de 
largeur nulle. 


Les rectangles R£ doivent être pris en compte dans la construction, même si leur aire est nulle, car 
ils servent pour pouvoir poser dessus à gauche un rectangle d’aire non nulle à l’étape suivante. 


Les rectangles sont indicés dans l’ordre décroissant de leurs aires, et non dans l’ordre de leur 
apparition dans la construction. 


2 Cet algorithme a été découvert par le mathématicien anglais William BROUNCKER (1620-1684) en 1668. 
https://fr.wikipedia.org/wiki/William_Brouncker 
https://www.pedagogie.ac-nantes.fr/mathematiques/enseignement/analyse-autour-de-In-1345111.kjsp 


5) Suites numériques : calcul de la valeur approchée de In(2) TermG spé maths © 2021 G.Lescaux 19 


La méthode de Bruncker peut être résumée par le schéma ci-dessous. 


À partir d’un rectangle R = ABCD dont le sommet supérieur droit C est sur l’hyperbole #4 
représentant la fonction inverse, on construit à l’étape suivante le rectangle R' = DMC'D' et le 
rectangle aplati R” = MCCM dont les sommets supérieurs droits respectifs C” et C sont sur 
l’hyperbole #, et dont le sommet M est le milieu de [DC]. 


Remarque : le rectangle R — ABCD peut être éventuellement aplati. Dans ce cas A = D, B=C, 
R = DCCD. 


ordonnées 


cf F. ) 
2 x+x 


% Xa + X X) abscisses 


Appliquons cette méthode à partir du rectangle R, tel que x, = 1, x2 = 2, y = 0. 
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n étape rectangle R}y 
étape longueur =| largeur = aire Ày = Ly x Ly 
base L, |hauteur #4 
1 On construit R, qui repose | L, = 1 FES 1 res 11 1 
sur l’axe des abscisses. T2 Hi 2 2. 1x2 
2 Sur R;, on pose R, et R£. 1 L LT 1? L 
LL; == Pr == A; =—X — = — = 
EE 
3 Sur R;, on pose R, et R3. 
i L: = — PS = — A3 = —X — = — — 
- - : 4 - 15 : 4 15 30 5 x 6 
Sur R;, on pose R, et Ra. = lp=— = 
; À LE f 14 26 7 x 8 
4 Sur R:, on pose R et Rs. 
RE 
| | : 8 : 48 : E 45 20 9 x 10 
Sur R;, on pose R, et R4. = be | A==x== == 
; : + : ”. Le 11 : 12 
Sur R,, on pose R; et R7. 
Le = — Pi = — Az = —X — = — — 
- . 7 8 ‘ 1 “ 8 1 182 13 x 14 
Sur R,, on pose Ra et Rs. 
Le == fo=— | === —— 
878 87 30 87 8° 30 240 15 x 16 
5 Sur Rs, on pose R, et RS. : m5 ? _ 8. A 1 
°7 16 | 7 153] ° 7 167 153 306 17 x 18 
Sur R£, on pose Ro et R6. ; CL ? m2 A 1 
07136! "7 95/7 16" 95 380 19 x 20 
Sur Re, on pose R:, et R°.. 1 8 L 8 1 1 
be lue Mar Ce — 
16 231 16 231 462 _21X 22 
Sur R£, on pose R:: et R9.. : LT r _ 2 A ROUES 1 
7716! 17 69) 127 16 69 552 23 x 24 
Sur R>, on pose R13 et Ro. : 7 r 2), AUS ot À 
177161! 1 325] 17 16 325 650 _25 x 26 
Sur RS, on pose R4 et R94.  — + . — +. 7 Le ee 
#77 16] 17 189/ 17 16 189 756 27 x 28 
Sur Re, on pose Ris et R9. ; _ 9 _"); +, 8, 4 1 
: : : 16 435 _ 16 435 870 29 x 30) 
Sur R3, on pose R.;ç et R;4. 
P 16 Lie pe me Re 
16 62 16 62 _992 31 X 32 
somme des 16 aires = 0,6777662022 
La somme de Brouncker de 16 premières aires vaut 
RL D LL LE 
+ 


380 462 552 650 756 870 992 
0.6777662022 
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Elle est égale à la somme de 2 qui vaut : 


S, (2) = EN AE ++ 
n (2) = An Vu + 3 n+n 


Avec ñn = 16, 


17 "18 19" 20 T1 va Ms 7 Tu TA 277 28 7 29 V30 "31 32 


0.6777662022 


L'égalité entre la somme de Brouncker et la somme de Riemann est vraie pour tout n € N°. 
Avec n = 10, la somme de Brouncker des 10 premières aires vaut 


Lo 2 2 1 1 Re 


2°12 30 56 90 132 182 240 306 380 
0.6687714032 


Elle est égale à la somme de Riemann, qui vaut : 


2 =) ——- PR tetes ts 
S10(2) Z10+k 11 12 13 14 15 16 17 18 19 20 


0.6687714032 
448372820160 


670442572800 
0.6687714032 


La fraction irréductible est : 


155685007 
232792560 
0.6687714032 
155685007 31 x 5022097 
232792560 2x 32X5X7X11 
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https://www.ma-calculatrice.fr/fraction-irreductibles-simplifier.ph 


Ma Calculatrice toute simple: Rendre une fraction irréductible 


448372820160 


670442572800 


RESET SIMPLIFIER SAUVEGARDE 


Résultat 


155685007 
2320792560 


Le nombre de rectangles construits vaut : 
- àlétapel1:1 
- à l’étape 2 : 2° =1 
- à l’étape 3 : 21 =2 
- à l’étape 4 : 22 = 4 
- à l’étape 5 : 2% =8 
- à l’étape 6 : 2* = 16 


- à l'étape m: 272 


Pour tout m € N°, le nombre total de rectangles construits en m étapes (c’est-à-dire avec m — 1 
étapes supplémentaires après l’étape 1) vaut : 


m—2 1— pm-1 
N= TN == 1 (2m) TL 42 = 2m 
k=0 
N = 281 


InN =In(2%-1) = (m —1)In2 


_nmN 


—1-— 
dé In 2 


— 2m-1 


Pour tout m € N°, le nombre total d’étapes pour construire N rectangles vaut : 
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J ESS 16 8 8 16/17 


8 l tu Le 1 1: 
at rs m0. 807 9T 153 


8/9 
8 | 4 EL BIT 
à SE 16231 462 16/19 
LE RE ET 475 
4 2 NE RARE ee 
GE efars eo | sen rss 16/23 
En ee 1 8 Î 2/3 
mi D 16*335 670 16/25 
ti ar É à 9 
see L Ent 189 4/7 
0.4 
0.3 
02 
0.1 
CPL ENNES ÉDAS EUR EZT 7 ES ENS EST 
4.6 8|16 2 16|8 16 4 6 ali6 
0.8 09 . 11 112 113 14 15 1.6 17 118 19 2 
A 2 = 1 ï 
-0.1 


Figure : construction des rectangles selon la méthode de Brouncker 


On a écrit en noir les coordonnées des points, 


en bleu les longueurs, en vert les largeurs, en violet les aires. 
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17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 


17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 
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0.9 


0.8 


0.7 


0.6 


0.5 


0.4 


0.3 


0.2 


0.1 


18 15/16 17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 x 116 


18 15/16 17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 


étape 1 
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17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 


étape 2 


0.3 


0.2 


0.1 


0 1—— 
18 15/16 17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 ges 116 


étape 3 
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0.9 
4/5 L 
—. 5, à, à À 
8 é° 1932. Ailrk 12 
6 9/11 
#7 Re 
0.6 
0.5 
0.4 
0.3 Li - 1 
2 1*x2 

0.2 
0.1 
0 À B, 
18 15/16 17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 TE 

étape 4 


9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 


étape 4 
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0:9 


07 


0:6 


1 
15 x 16 


C,1/2 


718 15/16 17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 7/4 29/16 15/8 31/16 
étape 4 
y 
+, 8. ‘À _  i 
\ 167 153 306 17x18 
Do Co 16/17 
La somme des aires des 16 restangles vaut 0.677766202207527 

0.9 


sie de à 
95 380 19 x 20 
D Ci016/19 


1 1 
90 9%x10 
0.8 


0.7 


1 1 


650 — 25 x 26 


0.6 Le Cu 16/27 15 = = NX 

Su D car 1,1 1] 3 
Les C5 16/29 ‘ 
D TRS #5 
1 Bis p, Ste 10/1 
15x16 16] [nc 


0.5 
B4 B3 BP T 
718 15/16 17/16 9/8 19/16 5/4 21/16 11/8 23/16 3/2 25/16 13/8 27/16 714 29/16 15/8 31/16 p 33/16 


étape 5 
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d) Algorithme du calcul approché de In(2) par le produit de 


Seidel, les sommes de Riemann et la somme de Brouncker 


Pour tout n € N”, la valeur de In 2 & 0,6931471806 peut être approchée par : 


e le produit de Seidel 
n 


_ 2 2 2 2 2 
MD = | = 


k=1 
= X 2 X - X X - 
RTE ee X E7rr 
DCE Vv2 
e les sommes de Riemann 
n-1 
D=y EE 
Tr An+tk on n+1 n+2 n+3 n+n-1 
ñn 
Or 
°n An+tk n+1 n+2 n+3 n+n 
e la somme de Brouncker 
ñn 
= ( 1 =: 1.1 1,1 1, . 1 1 
F7 AUk-1 2k/ 2 3 4 5 6 2n—1 2n 
> 1 D 1 
 AGk-Dx2k 1x2 3x4 5X6 (2n—1)x2n 


L’approximation est d’autant plus précise que n est grand. 
lim CL) = lim (2)= lim s,(2) = lim b,(1) = In2 = 0,6931471806 
n—+00 n—+00 n—+00 n—+00 


La différence entre s, (2) et b,, (1) est nulle pour tout n € N* : 


b\(1) — s,(2) = À ER ” —— 


: : (n +k) (2k—1) x 2k 
Des a 


k=1 

CCn+k-Gk-1)x2k] | n+k-—4k2+2k 
=> f-nxaein Dee 0) 
+ n + 3k — 4k? 
re 
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1. Écrire en langage universel, puis en langage Python, un algorithme In 2 (n) qui, pour tout 
n € N°: 


renvoie les quatre valeurs u, (1), r, (2), s, (2), b, (1) ; 

calcule leurs différences respectives à leur limite commune égale à In? ; 

donne la valeur plus précise de In 2 ; 

calcule la différence b, (1) — s,, (2), qui doit être nulle. On calculera cette différence de 
deux façons différentes : soit en calculant directement b,, (1) — s, (2), soit en calculant 
cette différence comme la somme de n termes : 


. n + 3k — 4k2 
; x — 1) x 2k(n +k) 


k=1 


2. Exécuter cet algorithme In 2 (n) pour n = 2* = 16. Parmi les quatre valeurs u,, (1), r, (2), 
Sn(2), ba(1), laquelle est la plus proche de In 2 & 0,6931471806 ? 
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Solution : 1. Algorithme en langage universel 


Fonction In 2 (n) 
=?» 
0—r 
Des 
0—b 
0O<— d 
Pour k allant de 1 à n Faire 


2 
PPX TR 
FinPour 


Afficher : « Le produit de Seidel P,, qui est le produit de —_— pour X allant de 1à n, 
vaut », p 
Afficher :«m2—P,=>»,In2-p 


Pour k allant de 0 à n — 1 Faire 


“A U 
T=T —— 
er 


FinPour 

Afficher : « La somme de Riemann R,, qui est la somme de — pour k allant de 0 à n — 
1, vaut »,r 

Afficher :««m2—R,=>»,Inm2-7r 

Pour k allant de 1 à n Faire 


L 
= S — 
$ Po 


FinPour 

Afficher : « La somme de Riemann S,, qui est la somme de =. pour k allant de 1 à n, 
vaut », S 

Afficher : «n2—S, =>»,In2—s 

Pour k allant de 1 à n Faire 


b = b+— 


(2k-1)x2k 
FinPour 


Afficher : « La somme de Bruncker B,, qui est la somme de pour k allant de 1 à 


es 
(2k—-1)X2k 
n, vaut », b 

Afficher : «In2 —B, =»,In2 —-b 


Pour k allant de 1 à n Faire 


n+3k-4k2 
d=d+ (2k—1)xX2k(n+k) 


FinPour 
Afficher : « La différence B, — S, est égale à », d, « ou », b —s 
Afficher : « La valeur plus précise de In 2 est », In 2 
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Algorithme en langage Python 


# Créé par user, Le 10/08/2021 avec EduPython 
from lycee import * 


def In2(n): 
p=1 


Q ou 
WOOU N l 
© © © © 


for k in range(1,n+1): # L'entier k varie de 1 à n 

p=p*2/(1+2**(1/(2**k))) 
print('produit de Seidel = P_‘",n,'= produit de 2 / [1 + 2**(1/2**k)] pour k = 1 à',n,'vaut',p) 
print('in(2) - P_',n,'=",ln(2)-p) 


for k in range(n): # L'entier k varie de 0 à n-1 

r=r+1/(n+k) 
print('somme de Riemann R_',n,'= somme de 1/(k+',n,') pour k = @ à',n-1,'vaut',r) 
print('in(2) - R_',n,'=",ln(2)-r) 


for k in range(1,n+1): # L'entier k varie de 1àn 

s=s+1/(n+k) 
print('somme de Riemann S_",n,'= somme de 1/(k+',n,') pour k = 1 à',n,'vaut'},s) 
print('lin(2) - S_‘,n,'=",ln(2)-s) 


for k in range(1,n+1): # L'entier k varie d_ 1àn 

b=b+1/((2*k-1)*2*k) 
print('somme de Brouncker = B_",n,'= somme de 1/[(2k-1)*2k] pour k = 1 à',n,l'vaut',b) 
print(‘'in(2) - B_',n,'=",ln(2)-b) 


for k in range(1,n+1): # L'entier k varie d_ 1àn 
d=d+(n+3*k-4*k**#2)/((2*k-1)*2*k*(n+k)) 

print('La différence B_",n,'moins S_",n,'est égale à',d,'ou',b-s) 

# La différence entre Les deux sommes précédentes est B n - S n = 0 


print('La valeur plus précise de 1In(2) est',ln(2)) 
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Algorithme en langage Python 


# Créé par user, le 10/08/2021 avec EduPython 


from lycee import * 


def In2(n): 
p=i 
1=0 
s=0 
b=0 
d=0 


for k in range(1,n+1): # l'entier k varie de 1 àn 

p=p*2/(1+2**(1/(2**k))) 
print('produit de Seidel = P_',n,'= produit de 2 / [1 + 2**(1/2*#k)] pour k = 1 à',n,'vaut',p) 
print(In(2) - P_',n,'=',In(2)-p) 


for k in range(n): # l'entier k varie de O0 à n-1 

r=1+1/(n+k) 
print(somme de Riemann R_',n,'= somme de 1/(k+',n,') pour k = 0 à',n-1,'vaut';r) 
print(In(2) - R_',n,'=',In(2)-r) 


for k in range(1,n+1): # l'entier k varie de 1 àn 

s=s+1/(n+Kk) 
print('somme de Riemann $S_',n,'= somme de 1/(k+',n,') pour k = 1 à',n,'vaut',s) 
print(In(2) - S_',n,'=',In(2)-s) 


for k in range(1,n+1): # l'entier k varie de 1 àn 

b=b+1/((2*k-1)*2*k) 
print('somme de Brouncker = B_',n,'= somme de 1/[(2k-1)*2k] pour k = 1 à',n,'vaut',b) 
print(In(2) - B_',n,'=',In(2)-b) 


for k in range(1,n+1): # l'entier k varie de 1 àn 
d=d+(n+3*k-4*k**2)/((2*k-1)#2#k*(n+k)) 

print('La différence B_',n,'moins S_',n,'est égale à',d,'ou',b-s) 

# La différence entre les deux sommes précédentes est B_n - S_n =0 


print('La valeur plus précise de In(2) est',In(2)) 
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2. Exécution de l’algorithme In 2 (n) pour n = 2* = 16 


>>> 1n2(2**4) 

produit de Seidel = P__ 16 = produit de 2 / [1 + 2**(1/2**k)] pour k = 1 à 16 vaut 0.6931508461384654 
in(2) - P__ 16 = -3.6655785201622493e-06 

somme de Riemann R_ 16 = somme de 1/(k+ 16 ) pour k = @ à 15 vaut 0.7090162022075267 

Iin(2) - R__16 = -0.015869021647581416 

somme de Riemann S_ 16 = somme de 1/(k+ 16 ) pour k = 1 à 16 vaut @.6777662022075267 

in(2) - S_ 16 = 0.015380978352418584 

somme de Brouncker = B_ 16 = somme de 1/[(2k-1)*2k] pour k = 1 à 16 vaut 0.6777662022075267 
in(2) - B__16 = 0.015380978352418584 

La différence B_ 16 moins S_ 16 est égale à 4.163336342344337/7e-17 ou 0.@ 

La valeur plus précise de In(2) est @.6931471805599453 


Sin = 2* = 16, b,(1) représente la somme des aires des 2° = 16 rectangles de R; à R46, qui 
donne une approximation de In 2. 


Ces rectangles sont construits en 4 étapes supplémentaires à partir du rectangle initial R:. 


1, 1, 1, 1, 1 1 1 EE 


Le EE LL D, 
*380 462 552 650 756 870 992 


0.6777662022 


In 2 
0.6931471806 


Parmi les quatre valeurs u,(1), r, (2), s,(2), b,(1), c’est le produit de Seidel u,(1) qui est le 
plus proche de In 2 & 0,6931471806. 
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e) Algorithme: fonction récursive utilisant la méthode de 
Bruncker 


La méthode de Bruncker peut être résumée par le schéma ci-dessous. 


À partir d’un rectangle R = ABCD dont le sommet supérieur droit C est sur l’hyperbole H 
représentant la fonction inverse, on construit à l’étape suivante le rectangle R' = DMC'D' et le 
rectangle aplati R" = MCCM dont les sommets supérieurs droits respectifs C”’ et C sont sur 
l’hyperbole #, et dont le sommet M est le milieu de [DC]. 


Remarque : le rectangle R — ABCD peut être éventuellement aplati. Dans ce cas A=D,B=C, 
R = DCCD. 


ordonnées 


cf, F ) 
2 "Xi + X 


x) Xi +2 X) abscisses 
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On considère un rectangle R = ABCD, dont les deux sommets inférieurs sont les points de 
coordonnées A(x,; y) et B(x:; y), et dont les deux sommets supérieurs sont les points de 


Æ 1 1 
coordonnées C (es =) et D Ga | =). 
2 


X2 
On considère la fonction récursive brouncker(x;, x, y, c) suivante, écrite en langage Python, qui 
renvoie, la somme de l’aire du rectangle R et de l’aire de tous les rectangles construits au-dessus 
de lui en c étapes supplémentaires, où c € N. 


Algorithme en langage Python 


def brouncker(x1,x2,y,c): 
On part d'un rectangle R ayant 
pour sommets inférieurs Les points de coordonnées (x1,y) et (x2,y) et 
pour sommets supérieurs Les points de coordonnées (x1,1/x2) et (x2,1/x2). 
Son aîre vaut Le produit de sa base (x2-x1) par sa hauteur (1/x2-y). 
On veut renvoyer La somme de son aire et de L'aire de tous Les rectangles 
construits au-dessus de Lui en c étapes supplémentaires. 
aire=(x2-x1)*(1/x2-y) 
Af c>0: 
aire=aire+brouncker(x1,(x1+x2)/2,1/x2,c-1) 
on ajoute L'aire du rectangle ayant 
pour sommets inférieurs Les points de coordonnées (x1,1/x2) et ((x1+x2)/2,1/x2) et 
pour sommets supérieurs Les points de coordonnées (x1,2/(x1+x2)) et ((x1+x2)/2,,2/(x1+x2). 
Son aire vaut Le produit de sa base (x1+x2)/2-x1 par sa hauteur 2/(x1+x2)-1/x2. 
aire=aire+brouncker((x1+x2)/2,x2,1/x2,c-1) 
on ajoute L'aire nulle du rectangle aplati ayant 
pour sommets inférieurs Les points de coordonnées ((x1+x2)/2,1/x2) et (x2,1/x2) et 
pour sommets supérieurs Les points de coordonnées ((x1+x2)/2,1/x2) et (x2,1/x2). 
Son aire est nulle : elle vaut Le produit de sa base x2-(x1+x2)/2 par sa hauteur 1/x2-1/x2 = 
return aire 
on retourne La somme de L'aire du rectangle R initial et de L'aire de tous Les] rectangles 
construits au-dessus de Lui en c étapes supplémentaires. 
Si on part du rectangle initial x1=1, x2=0, y=0, 
Le nombre total de rectangles obtenus avec c étapes supplémentaires vaut n = 2**c 
brouncker(1,2,0,c) = B_n renvoyé par lLn2(n) avec n = 2**c donc c = lLn(n)/tln(2) 


+ He # He # # 


# # # # # # # # 


# H # # # 


def brouncker(x1,x2,y,c): 
# On part d'un rectangle R ayant 
# pour sommets inférieurs les points de coordonnées (x1,y) et (x2,y) et 
# pour sommets supérieurs les points de coordonnées (x1,1/x2) et (x2,1/x2). 
# Son aire vaut le produit de sa base (x2-x1) par sa hauteur (1/x2-y). 
# On veut renvoyer la somme de son aire et de l'aire de tous les rectangles 
# construits au-dessus de lui en c étapes supplémentaires. 
aire=(x2-x1)*(1/x2-y) 
if c>0: 
aire=aire+brouncker(x1,(x1+x2)/2,1/x2,c-1) 
# on ajoute l'aire du rectangle ayant 
# pour sommets inférieurs les points de coordonnées (x1,1/x2) et ((x1+x2)/2,1/x2) et 
# pour sommets supérieurs les points de coordonnées (x1,2/(x1+x2)) et ((x1+x2)/2,2/(x1+x2). 
# Son aire vaut le produit de sa base (x1+x2)/2-x1 par sa hauteur 2/(x1+x2)-1/x2. 
aire=aire+brouncker((x1+x2)/2,x2,1/x2,c-1) 
# on ajoute l'aire nulle du rectangle aplati ayant 
# pour sommets inférieurs les points de coordonnées ((x1+x2)/2,1/x2) et (x2,1/x2) et 
# pour sommets supérieurs les points de coordonnées ((x1+x2)/2,1/x2) et (x2,1/x2). 
# Son aire est nulle : elle vaut le produit de sa base x2-(x1+x2)/2 par sa hauteur 1/x2-1/x2 = 0 
return aire 
# on retourne la somme de l'aire du rectangle R initial et de l'aire de tous les rectangles 
# construits au-dessus de lui en c étapes supplémentaires. 
# Si on part du rectangle initial x1=1, x2=0, y=0, 
# le nombre total de rectangles obtenus avec c étapes supplémentaires vaut n = 2**c 
# brouncker(1,2,0,c) = B_n renvoyé par In2(n) avec n = 2**c donc c = In(n}/In(2) 
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Questions 


1. Quel résultat renvoie l’exécution de brouncker (2 ; 3 =, =. , 2) 7 
Détailler la procédure suivie par l’ordinateur pour bte ce résultat. 

Que représente ce résultat ? 

2. Pour tout c € N, que représente le résultat renvoyé en exécutant brouncker(1,2,0, c) ? 
Justifier que pour tout c € N, sin = 2°, l’algorithme In 2 (n) vu dans la partie précédente 
donne pour b, (1) le même résultat que brouncker(1,2,0, c). 

3. Que vaut brouncker(1,2,0, c) si c = 4 ? Que représente ce résultat ? 

Vérifier que si c = 4 donc n = 2* = 16, l’algorithme In 2 (16) donne pour b46(1) le 
même résultat que brouncker(1,2,0,4). 


Solution 


1. Exécutons brouncker G, . ,2) 

| 7 3 4 1 1 
cre-(f-De 6 
2>0 


3 4 3 13 à 
On ajoute à aire la valeur brouncker Ë 17 1) = brouncker G, ai 1) 


3,7 
On ajoute à aire la valeur brouncker (24. 1) = brouncker (£, 2,2,1) 
2 4 7 8 4 7 
13 4 
Exécutons brouncker , ,1) 
2° 8 ty 
| (+ =) x É =) = 1 . 
FO AD à 1 D Me 
Le0 
3,13 
On ajoute à aire la valeur brouncker (: 40) brouncker , _ ,0) 
2 13 2°16/13° 


| = (© )x (Se =) = 1 — 
FO M 2) 02 a 60 
Exécutons brouncker (2 À : , 1) 


| = ( x G =) = 0 
MENT à 7 7) 


1>0 


1957 

On ajoute à aire la valeur brouncker (2. 520) - brouncker (Æ, 22,0) 
| É =) œ -) À . 
=———|K———| = — — 
Fe Né 6) Ur d'en. À 
brouncker (£, . = 2) renvoie 
24/2 
Aa + A7 + As + Aa 


C’est la somme de l’aire du rectangle R, initial et de l’aire des trois rectangles R;, R;3, R:4 
construits au-dessus de lui en 2 étapes supplémentaires. 
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>>> brouncker(3/2,7/4,1/2,2) 
0.026212861212861212 

>>> 1/56+1/182+1/650+1/756 
0.026212861212861212 


2. Pour tout c € N, l’algorithme brouncker(1,2,0, c) renvoie, à partir du rectangle R, de base 
1 1 4 
2—1=1et de hauteur — 0 = 7 la somme de son aire ; © de l’aire de tous les rectangles 


construits au-dessus de lui en c étapes supplémentaires. Avec un total c + 1 étapes, on obtient un 
total de n = 2° rectangles dont la somme des aires vaut brouncker(1,2,0, c). 


Donc, pour tout c € N, si n = 2°, l’algorithme In 2 (n) donne pour b, (1) le même résultat que 
brouncker(1,2,0, c). 


Ce résultat donne une approximation de In 2, qui est d’autant plus précise que c est grand (donc 
n est grand). 


3. Si le nombre d’étapes supplémentaires vaut c = 4, on obtient 


>>> brouncker(1,2,@,4) 
0.6777662022075268 


C’est la somme des aires des 2* = 16 rectangles de R; à R:6, construits en 4 étapes 
supplémentaires à partir de R,. Cette somme donne une approximation de In 2. 


1 1 1,1, 1 1 1 1 1 
212 30 56 90 132 182 240 306 


*380 462 552 650 756 870 992 


0.6777662022 


In 2 
0.6931471806 


Si c = 4 donc n = 2* = 16, l’algorithme In 2 (16) donne pour b46(1) le même résultat que 
brouncker(1,2,0,4). 


>>> 1n2(2**4) 

produit de Seidel = P__ 16 = produit de 2 / [1 + 2**(1/2**K)] pour k = 1 à 16 vaut 0@.6931508461384654 
in(2) - P__16 = -3.6655785201622493e-06 

somme de Riemann R__ 16 = somme de 1/(k+ 16 ) pour k = @ à 15 vaut 0@.7090162022075267 

in(2) - R__16 = -0.015869021647581416 

somme de Riemann S_ 16 = somme de 1/(k+ 16 ) pour k = 1 à 16 vaut @.6777662022075267 

in(2) - S__16 = 0.015380978352418584 

somme de Brouncker = B_ 16 = somme de 1/[(2k-1)*2k] pour k = 1 à 16 vaut 0.6777662022075267 
in(2) - B__ 16 = 0.015380978352418584 

La différence B_ 16 moins S_ 16 est égale à 4.163336342344337e-17 ou 0.0 

La valeur plus précise de In(2) est 0.6931471805599453 


>>> brouncker(1,2,0,4) 
0.6//7/7662022075268 
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f) Annexe : égalité de la somme de Brouncker et de la somme de 


Riemann 
Théorème 
Pour tout n € N°, on considère la somme de Brouncker b,, et la somme de Riemann s,, définies 
par : 
DT -n)d 
0 ÏJ 2k—1 2k/ Z(2k-1)x2k 
j=1 k=1 k=1 
et 


2n n 
#2: ee 
i + k 
i=n+ k=1 
Alors on a les propriétés suivantes : 


Dh=s=; 


2) Pour tout n € N*, 
1 1 1 


2n+1 2n+2 ” (2n+D(G@n+2) 


bin se 


3) Pour tout n € N”, 
b;=S, 


Démonstration 
Vérifions d’abord l’égalité des différentes écritures de b, et de s,.. 


oo sn  : 1 
Te E 2 504 "5 6 2n-1 ?2n 


“Es: 


L 1 1 1 


5 a 


Une autre écriture de cette somme est 


En effet, 
2n-1 : 2n— à 2n— ; 2n—2 : 2n— ; 
MN CDN ED STE SE N CD 
FT IL Gil D i+1 Ai+i ln ET 
i=0 i=0 i=0 i=0 i=1 
i pair i impair i pair i impair 
: (—1)2k-2 - (—1)2#-1 : . (—1)2* : (—1)2* 
Like 24i 2k—1+1 A2k-1 2k 
k=1 k=1 k=i k=1 
ñn n n ñn 
Dr Da x) 
| Zi2k—1 2k 2k—1 2k/ Ai(2k—-1)X2k 
k=1 k=1 k=1 k=1 
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La somme de Riemann vaut 


> -  — LA à 
: _ is 
" Antk n+1 n+2 n+3 n+n 


En effectuant le changement de variable i = k + n dans s,,, k varie de 1 à n, donc i = k + n varie 
de1+nàn+n = 2n. Donc une autre écriture possible de cette somme est 


1) Montrons que b; = S1 = : 


2x1 ' 

1 J+1 | 1+1 = 2+1 1 _—1 1 1 
— J 
J=1 


1 2 1 2 2 2 
2x1 2 
: EN 
AT mn 2 
i=1+1 i=2 
b, = _1 
RS 
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2) Montrons que pour tout n € N’, 


1 1 1 


b — bb, =Ss —S D DA Çn + 1)(2n +2) 
n+1 n n+1 n Ton + 1. 2n +2 + 2 |” (2n + 1)(2n + 2) 


Preuve : 
n 
1 1 
ED Cr 2) 
2k—1 2k 
k=1 
n+1 n 
PE n° a a an 
HT 2k=-1. Li 2k=1. 2n+1)-1 2(n+1) 
2. 1 n. 1 1 
Tondrel | 2n+2 2n+1 2n+2 
: ; 1 1 1 
"HR 2n+1 2n+2 (2n+1)(2n+2 
ou bien 
2n (—1)/#1 
— 
j=1 J 
ASTOUS _ (= D C= D" (— DRE ( LIFATÈTE 
CD DE eme un 
2n +1 2n +2 
ee È Dsl 1 1 
= D, ++ —— > pb, — 
2n +1 2n +2 2n+1 2n+2 
1 1 1 
D1 b, = 


2n+1 2n+2 ” (n+1D(@n+2) 
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n+1 n+1 


1 
nn) n+k+1 
k=1 k=1 


On effectue le changement de variable i = k + 1. Lorsque k varie de 1 à n + 1, i = k + 1 varie 
de1+1=2àn+1+1-n+2. 


n+1 n+2 n+2 n+2 


1 1 1 " 1 . 1 1 : 1 
S. = — = ———© = — — —— —— ——— 
Ee D » n+i n+1 n+1 » n+i n +1 > n+i 
k=1 i=2 i=2 i=1 


n 


1 1 1 1 
n +1 in+i n+n+i n+tn+2 


1 PRE. dd 2x1 
CET EL 2n62 nel © '2nil 242 2EMmID 
1= 
,_17 1 , 17 1 2 : 1 , 172 
= S. — = $ 
M 2n+1 2n+2 2(n+1) " 2n+1 2n+2 
1 1 
n + — 
2n +1 | 2n+2 
1 1 | 
S —$,, = —— = = ————— 
MAUR 2n+1 2n+2 (2n+1)(2n+2) 
ou bien 
an 
1 
“à 
i=n+1 
2(n+1) 2n+2 2n+ 2n+2 
: » L _ Y == 1 DE + 
SR D Zi mt ‘0 PEN 
i=n+1+1 i=n+2 à i=n+ i=n+1 


1 1 1 1 
Ge a nTon+1 2n+2 


i=n+ 
_ A a À à 2 
TATDREL Are 20 © 2nel 242 nil 
1 1 1 1 


À —— #25, + —_— — — _— 
2n+1 2n+2 * 2n+1 2n+2 


1 1 1 
S — S, = ————— =] —_. — …—] . _— 
TOR 2n+1 2n+2 (2n+1)(2n+2) 


Conclusion : 
1 1 1 


b... — b. = _ a — OnilOni2) 
n+1 Dr Snt1 Sn 311 2n+2 ” Cn+1(2n+2) 
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3) Montrons que pour tout n € N’, 


b, — Sn 
On a montré que : 
1 
1) b = s1 = 2 
2) Pour tout n € N*, 
1 1 1 


babes es] — ——— ———— 
PNR CT SE EL 22 CrEFUORET 
Il en résulte que pour tout n € N”, 
Dn+1 — Snt1 = Dn — Sn = bi — S1 = 0 
Donc 
bn = Sn 
Voici une démonstration directe du point 3 (sans utiliser les points 1 et 2) 


http://ww2.ac-poitiers.fr/math/sites/math/IMG/pdf/H-Tarfaoui-9.pdf 


Pour tout n € N°, 


: 1 3 i 11,11 1 1 
 AUk-1 2k 2 3 4 5 6 2n—-1 2n 
_ 1 1,1,1 | 1,1,1 1 1,1. - 1 té + À 
E D 2 1 5 #4 4 4.5 6 © 6 2n—-1 2n 2n ?2n 
” 1 1,1,1 1 1,1,1 1 1,1, L 1 TL À  : 
_ à 2 9 4 À. à #4 5 6 6 € 2n—-1 2n 2n ?2n 
ai ii + : { ,1 4 ET À Æ À Lt. + 
_ 3 4 5 6 2n—-1 2n 2 2 4 4 6 6 2n 2n 
=. du + - Le 2(S + +24 + 
_ 2 3 à 5 6 2n—-1 ?2n 2 4 6 2n 
+ : i— (i+s+st +2 
_ 2 3 4 5 6 2n-1 ?2n + à n 
Sp pr en +. + + ++ Là (1+5+3+ +) 
DRE re n n+1 n+2 n+3 2n—1 ?2n 2°: n 
= 1 . TL T (: 1. À =) 
A D = Di D: 
1 1 1 1 1 + 1 


ntitn+2 ns 2n-1 2m Link 
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Application de l’égalité b, = s, 

Exercice 1 

XXIe Olympiades internationales de mathématiques : exercice 1 du sujet du 2 juillet 1979 
http://imo-official.com/problems.aspx 


Soient a et b des entiers strictement positifs vérifiant : 
a 1 1 1 1 1 1 1 1 


bo 2 3 4 5 0e ‘ou ci vo 
Montrer que 1979 divise a. 
Solution 


Pour tout n € N*, 


bis. 
où 
n 
= Ù( 1 ) = 1 1 en 1 1 1 
* AVk-1 2k/ 2 3 4 5 6 2n—1 2n 
et 
n 
=). A ++ 
7% Lin+k n+1 n+2 n+3 n+n 
En particulier, si n = — = 659, 
béso = S659 


: 1 
La somme de Brouncker augmentée de . 


de a Le 
#7 1319 ZVk-1 2k/ 2 3 45 6 1317 1318 1319 D 


pa N . » 1 
est égale à la somme de Riemann augmentée de . 


659 


à 
L 1319 23659+k 660 661 662 1317 1318 1319 


- 4 . 4 : ++ 1 + - de - 
_1 660 661 662 1317 1318 1319 
9 1 1 
D _ — 
1319 1318 1317 662 661 660 
1 n 1 A 1 . i : 1 A 1 
_1 660 661 662 1317 T1 150 
- + + + did - + - + - 

1319 1318 1317 662 661 660 
= 
7 2\660 1319 661 1318 662 1317 1317 662 1318 661 1319 660 

1 1979 1979 1979 1979 1979 
= 
2 \660 x 1319 661 * 1318 662 x 1317 1317 X 662 1318 x 661 
p 


a | = 1079 %X— 
. 1319 % — q 
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où p et q sont des entiers naturels vérifiant 


[TS 


( 1 1 1 1 1 


= ———_—_—_—_— à + ——— 
q 2\660x1319 661xX1318 662 x 1317 1317 X 662 1318 x 661 


1 
F LL ) 


L'égalité biso = Sés9 implique l’égalité 

a P 

—= 1970 X — 

b q 
Donc 

axqg=1979xXpxb 

Donc 1979 divise a X q: 

ax q 

=pXxb 
1979 P 


Or 1979 est un nombre premier 
https://www.nombres-premiers.fr/liste.html 


https://www.nombres-premiers.fr/1979.html 


1979 ne divise pas q, car q est formé d’entiers naturels strictement inférieurs à 1979. 
1979 et q sont premiers entre eux 


On utilise le théorème de Gauss : 


Théorème de Gauss 
Soit à, b, c trois entiers relatifs, avec a et b non nuls. 
Si a divise le produit bc, et a et b sont premiers entre eux, alors a divise c. 


1979 divise le produit a x q, et 1979 et q sont premiers entre eux, donc 1979 divise a. 
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Exercice 2 (inspiré de l’exercice 1) 

Énoncé : http://ww2.ac-poitiers.fr/math/spip.php?article434 

Solution : http://ww2.ac-poitiers.fr/math/sites/math/IMG/pdf/H-Tarfaoui-9.pdf 
http://ww2.ac-poitiers.fr/math/sites/math/IMG/pdf/W-Mesnier.pdf 


Soient a et b des entiers strictement positifs vérifiant : 
a 1 1 1 1 1 1 1 1 


p=772t3 ts 6t"Tis20 1330! 1331 
Montrer que 1997 divise a. 
Solution 
Pour tout n € N°, 

bn = Sn 

où 

D, 1 1 À #1 à 1 1 

=) (=) Da br Ne, 

et 


1 1 1 
D a ea on 
En particulier, si n = = 665, 


bées = Sées 


u 1 
La somme de Brouncker augmentée de a 


best d (— —)=1 cr ru in Re 
1331 ZA Uk-1 2k/ 2 3 4 5 6 1329 1330 1331 D 
est égale à la somme de Riemann augmentée de . 

665 


EE +: 
$665 1331 4 665+k 666 667 668 1329 1330 1331 


=1 


RS et 
_1 666 667 668 1329 1330 1331 
Le + + ne 
1331 1330 1329 668 667 666 
EE _ — 
_1 666 667 668 1329 1330 1331 
a 
1331 1330 1329 668 667 666 
= RE  : es 
7 2\666 1331 667 1330 668 1329 1329 668 1330 667 1331 666 
1 1997 1997 1997 1997 1997 
= ss ee + 
2\666x 1331 667 x 1330 668 x 1329 1329 x 668 1330 x 667 
p 
———— )= 1997 x— 
a q 
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où p et q sont des entiers naturels vérifiant 


[TS 


( 1 L 1 1 1 


= ———_—_—— à + —— 
q 2\666x1331 667x 1330 668 x 1329 1329 x 668 1330 x 667 


£ 
u LL 4 =) 


L'égalité bis = Sécs implique l’égalité 

a P 

—— 1997 X— 

b q 
Donc 

axq=1997xpxb 

Donc 1997 divise a X q: 

axq 

=pXxb 
1997 P 


Or 1997 est un nombre premier 
https://www.nombres-premiers.fr/liste.html 


https://www.nombres-premiers.fr/1997.html 


1997 ne divise pas q, car q est formé d’entiers naturels strictement inférieurs à 1997. 
1997 et q sont premiers entre eux 


On utilise le théorème de Gauss : 


Théorème de Gauss 
Soit à, b, c trois entiers relatifs, avec a et b non nuls. 
Si a divise le produit bc, et a et b sont premiers entre eux, alors a divise c. 


1997 divise le produit a x q, et 1997 et q sont premiers entre eux, donc 1997 divise a. 
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